TITLE 
TYPE 
DURATION 
DESCRIPTION 



METHOD OF USE 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



Aux. 
ILLINOIS CODE M 7 - 150 



Principal Axes Factor Analysis (DOI only) 
Entire program, includes DOI 
Approximately the same as for M 3 - 117 

A modification of M 3 which prints the eigenvectors scaled 
"by the square root of the corresponding eigenvalue (i.e., 
latent root). The eigenvectors "become factors by this 
scaling. Each factor is printed as a column terminated by 
the character "N". The eigenvalues are punched as a 
separate column after the last factor. The user of this 
program has to supply only the problem tap which is 
described below. 

(1) Read program tape in until it stops* 

(2) Insert problem tape in reader and input the symmetrical 
matrix and two parameters. 

(3) Replace program tape in reader and continue input. 
The first part of the computing will begin before the 
program tape has been read completely. When this part 

of the computing is finished the rest of the program tape 
will be input, computation completed, and the results 
punched out in the manner described above. 



PREPARATION OF 
PROGRAM TAPE 
PROBLEM TAPE 



TAPES 



Obtain a copy of the tape in the computer library. 
The first part of this tape consists of the matrix to be 
factored. The lower left triangular half of the matrix is 
punched by rows, each element being in the form required 

by Illinois Code N 3, 23, Decimal Number Input. Each 

-2 
element should be scaled by 10 before being punched. The 

last element in the matrix should be followed by the 

character "N". The second part of the tape consists of two 

decimal numbers, each followed by a sexadecimal character. 

If (d) is the number of places desired in the printed results 



-2- 



COIMEHTS 



EXAMPLE 



and (p) is the order of the symmetric matrix to be factored 

then the second part of this tape can "be represented as 

(d) J (p) F. 

The capacity of the routine is a matrix of order 23. The 

-1 



factor loadings will he scaled by 10 

v -2 



when printed and 
the eigenvalues by 10 . Either unities or communalities 
can be used in the leading diagonal. 
Let .906 .^38 .072 

.438 .891 -o309 
.072 -.309 0650 
be the symmetric matrix which we are to factor, securing 
results to four decimal places. Then we will make a copy 
of the program tape in the library and make the following 
problem tapes 
spaces 
+00906 
+00^38 
+00891 
+00072 
-00309 
+OO65O 
N 

5 
J 
spaces 



J 



the lower half of the matrix 
punched by rows 



the printed results will have a 
scaling zero and four significant 
figures 



the matrix is of order three 



The results will be printed as 
+07^95 
+O8563 
-02891 

N 
-02526 
+03231 



-3- 

+03022 

N 
+05296 
-02308 
+06892 

N 



\> 




factor loadings 
scaled by 10 " 



eigenvalues scaled by 10 



-2 



Coded by J. 0. Neuhaus 
November 195^ 



Retyped: 9/24/59 



LOCATION 



ORDER 



JJ2CBS. 



PAfiE 1 M 7 



Decimal Order Input - Code XI 



00 100K 



Input Decimal Fractions - Code N2 



10 



11 



12 



13 



li* 



15 



00 10K 




40 170F 


*" 


50 L 




26 100F 


- 


4l F 


from 9 


81 4F 


from 6 


LO 28L 




32 6L 




L4 28L 




50 F 




?4 28L 




S5 F 




40 F 




26 2L 




42 7L 


from 3 


L5 F 




40 ( )F 


by 6 


LI 10L 




40 10L 




32 1L 




50 4F 




L5 OOf 


by 8 


74 4f 




S5 F 




10 IF 




L4 29L 




40 7F 




42 17L 




L5 4F 




L4 4f 




40 5F 




4l IF 




F5 17L 





Read in correlation matrix 



Read in parameters and convert 



Form dependent parameters 



LOCATION 
16 

17 
18 

19 
20 

21 
22 

23 

24 

25 
26 
27 
28 
29 
30 



ORDER 

42 20L 
00 IF 
L5 30L 
40 ( )F 

F5 IF 

40 IF 
LO 4F 
32 26L 

41 F 

41 ( ) F 

F5 F 
40 F 
LO 4f 
32 24L 
F5 20L 

42 20L 
22 20L 
L5 17L 
F4 4f 
42 17X 
22 15L 
L5 4f 
4o 13F 
24 999F 
00 F 

00 10F 
00 F 
00 170F 
80 F 
00 F 
24 ION 



NOTES 
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by 13, 25 



t 



by 16, 23 
from 24 



Generate unit matrix 






from 22 



— ' Generate unit matrix 



from 19 




10 x 2 



-39 



first storage location -Constants 



-1 



Enter parameters and matrix 



LOCATION 







ORDER 



NOTES 
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M 7 



Save parameters 



00 20K 

Eigenvalue Routine as used for Code 117, modified. 
00 152KJ 

Square Root - Code Rl 
00 161K 
00 F 
00 S3 
80 F 
00 Sh 
00 F 
00 S7 
26 21N 
00 10K 
L5 161F 
kO 3F 
26 999F 
00 F 

126 ION 
00 110KI 
Single Column Print - Code P£ 
00 10K 



kl F 
50 F 
L5 (170)F 
i+0 (170 )F 
F5 F 
^0 F 
Kk F 
00 20F 
Fk 1L 
^0 1L 



Plant print parameter 



from 6, by i+. 
by U 



— Store Eigenvalues consecutively 



LOCATION 



ORDER 



amm 



NOTES 
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M 7 



8 



10 



11 



12 



13 



14 



15 



16 



17 



18 



19 



20 



L5 
LO 
36 
41 
L7 
50 
26 
40 

19 
L4 
46 

F5 
42 
F5 
40 
LO 
36 
L5 
42 
92 
41 
41 
50 
79 
22 
50 
26 

L5 
L4 

42 



F 

162F 

1L 

F 

(170 )F 

7L 

152F 

(80 )F 

18F 

7L 

7L 

8L 

8L 

r 

F 

162F 

7L 

163F 

16L 

171F 

4f 

5F 

(80 )F 

( )F 

17L 

17L 

110F 

l6L 

162F 

16L 



F5 5F 
40 5F 



by 10, from 13 

to sq, rt. 

by 11 



J 



- Take square root of eigenvalues 



Take square root of eigenvalues 



from 27 

from 21, by 26 

by 14, 19, 24 



to S. C. P. 



Punch factor loadings 

("N" punched after last loading 
in each factor) 



LOCATION 



21 
22 

23 
24 

25 

26 

27 
28 

29 
30 
31 
32 
33 
34 



ORDER 



LO 


162F 


36 16L 


92 


770F 


92 


139F 


F3 


4f 


40 4f 


L4 16 3F 


42 16L 


19 18F 


L4 l6L 


46 16L 


L5 


4f 


LO 


162F 


32 


15L 


92 


171 


92 


171F 


41 


5F' 


L5 


(170 )F 


22 


30L 


50 


30L 


26 


110F 


F5 29L 


42 


29L 


F5 


5F 


40 


5F 


LO 


162F 


32 


29L 


OF 


F 


26 


ION 



NOTES 
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by 32, from 34 



to S. C. P. 



Punch factor loadings 

("N" punched after last loading 
in each factor) 



Punch spaces 



— Punch eigenvalues 



< 11 11 1 



End 

Enter output code 



